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ABSTRACT 

Aims. When selecting flux-limited cluster samples, the detection efficiency of X-ray instruments is not the same for 
centrally-peaked and flat objects, which introduces a bias in flux-limited cluster samples. We quantify this effect in the 
case of a well-known cluster sample, HIFLUGCS. 

Methods. We simulate a population of X-ray clusters with various surface-brightness profiles, and use the instrumental 
characteristics of the ROSAT All-Sky Survey (RASS) to select flux-limited samples similar to the HIFLUGCS sample 
and predict the expected bias. For comparison, we also estimate observationally the bias in the HIFLUGCS sample 
using XMM-Newton and ROSAT data. 

Results. We find that the selection of X-ray cluster samples is significantly biased (~ 29%) in favor of the peaked, 
Cool-Core (CC) objects, with respect to Non-Cool-Core (NCC) systems. Interestingly, we find that the bias affects the 
low-mass, nearby objects (groups, poor clusters) much more than the more luminous objects (i.e massive clusters). We 
also note a moderate increase of the bias for the more distant systems. 

Conclusions. Observationally, we propose to select the objects according to their flux in a well-defined physical range 
excluding the cores, 0.2rsoo — »"500, to get rid of the bias. From the fluxes in this range, we reject 13 clusters out of the 
64 in the HIFLUGCS sample, none of which appears to be NCC. As a result, we estimate that less than half (35-37%) 
of the galaxy clusters in the local Universe are strong CC. In the paradigm where the CC objects trace relaxed clusters 
as opposed to unrelaxed, merging objects, this implies that to the present day the majority of the objects are not in a 
relaxed state. From this result, we estimate a rate of heating events of ~ 1/3 Gyr" 1 per dark-matter halo. 

Key words. Galaxies: clusters: general - X-rays: galaxies: clusters - Galaxies: clusters: intracluster medium 



1. Introduction 

Galaxy clusters are the biggest gravitationally-bound struc- 
tures in the Universe. They are filled with a hot (kT ~ 
1 — 10 keV) plasma, the Intra-Cluster Medium (ICM), which 
has been heated to X-ray emitting temperatures by grav- 
itational colla pse. In the original "cooling-flow" scenario 
(Fabian 1994), all clusters after relaxation from a major 



merging event should evolve into the cooling-flow state, 
where the gas condensates in the central regions and then 
cools through radiative processes until it eventually forms 
stars. This paradigm was supported by observations of 
prominent surface-brightness peaks and temperature drops 
in the cores of clusters, which were supposed to be associ- 
ated with the central cooling flow. However, this picture was 
not confirmed by the latest generation of X-ray telescopes 
{XMM-Newton, Chandra) , which found no spectroscopic 
evidence for the presence of the cooling gas in the cen- 
tral regions of clusters predicted by the coo ling-flow model 



(Peterson et al. 2001 



2001). These results 



Kaastra et al 

lead to the revision of the classification of galaxy clusters 
which were thereafter categorized i nto "cool-core" (CC) 



and "non-cool-core" (NCC) objects (Molendi & Pizzolato 
20011. Since in the center of clusters the cooling time can be 
much shorter than the Hubble time, these results imply the 
existence of a heating mechanism which is responsible for 
quenching the cooling flow. Feedback from Active Galactic 



Nuclei (AGN) in the central galaxie s is the most probable 
sourc e of heating in the ICM (e.g., McNamara & Nulsen 



2007) 



Apart from the nature of the heating source, the fail- 
ure of the cooling-flow model also has repercussions on the 
formation scenario of galaxy clusters. Indeed, some numeri- 
cal studies of cluster mergers have fo und that cool core s are 
hardly disrupted by merging events ( Poole et al.|20 08 ) , and 
that the state of a cluster (CC or N GC) is deter mined onc e 
for all during the formation process ( McCarthy et al.|20 08). 
Recent observational works have shown that the popula- 

which 



tion of clusters is bimodal (Cavagnolo et al. 
could support this scenario because of the 



2009J), 
ack of inter- 



mediate objects. However, other studies found no ev idence 
for bimoda lity in the distribution of clusters (e.g., Pratt 



et al. 2010), and the recent identification of regions rem- 



miscent of cool cores in merging cl usters ("cool-core rem- 



nants", |Rossetti fc Molendi (2010)), probably associated 
with ancient cool cores disrupted by merging events, ar- 
gues against this idea. Thus, the question of the formation 
process of cool cores is still an open one. 

An important clue to the understanding of the forma- 
tion process of cool cores is the observed ratio between 
CC and NCC clusters in the local Universe, which depends 
strongly on the different formation scenarios. Indeed, in the 
old cooling-flow scenario it was believed that the majority 
of clusters (70-90%) in the local Universe had a cooling flow 
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( Peres et al. 1998 ) , which is expected if all clusters naturally 
evolve into the cooling-flow state. In this case, one might 
also expect a strong dependence of the CC fraction on the 
redshift. Conversely, if the state of a cluster is determined 
only by an initial entropy injection event as suggested by 
some simulations (McCarthy et al. 2004), a very different 
behavior can be expected. 

Observationally, to measure the fraction of CC objects 
it is crucial to use a sample which is as complete and as free 
of selection biases as possible. In this respect, the very dif- 
ferent surface-brightness profiles observed in CC and NCC 
objects might introduce a bias in the selection of X-ray 
flux-limited samples. Indeed, CC systems exhibit a promi- 
nent surface-brightness peak, as opposed to NCC clusters 
which show flat emission profiles. The diversity of surface- 
brightness profiles might introduce a bias in favor of CCs in 
X-ray flux-limited samples. In this paper, we investigate the 
impact of such a flux selection on the measurement of the 
CC over NCC ratio when computed using a well-studied, 
nearby cluster sample, HIFLUGCS, which was extracted 
from the ROSA T All-Sky Survey (RASS) 

HIFLUGCS ( |Reiprich fc Bohringer|2002"t is a complete 
flux-limited sample extracted from RASS data. It comprises 
64 clusters with a flux higher than 2.0 x 10 -11 ergs cm~ 2 
s _1 (0.1-2.4 keV band). According to the authors, it is to 
the present day the largest complete X-ray selected sam- 
ple of galaxy clusters. Additionally, an extended sample of 
106 objects also exists, which comprises objects which were 
originally excluded for completeness reasons (in particular 
because of high absorption). The objects of the extended 
sample are nearby (z < 0.2) and span a luminosity range 
between ~ 10 43 and ~ 10 45 h^ 2 ergs s _1 . 

From ROSAT and ASCA data, 



Chen et al 



(2007) es 



timated the fraction of CC to be 49% in the HIFLUGCS 
sample. More recently, studies have been carried out us- 



ing Chandra data ( |Mittal et al.||2009| |Hudson et al. 12010 1 
which clearly demonstrate the existence of a class of inter 



mediate objects. In particular, Hudson et al 
after, H10) claim a fraction of 44% 
and 28% "weak CC" clusters, where the classification 



, P610| ) (here- 
"strong CC" objects 
of 



the objects was performed using the central cooling time 
(CCT), which was found to be the best indicator of the 
state of a cluster. Strong CC objects were defined as ob- 
jects which exhibit a central cooling time < 1 Gyr, while 
weak CC objects have 1 Gyr < CCT < 7.7 Gyr, and clusters 
with a CCT > 7.7 Gyr are classified as NCC. In addition, 
H10 also analyzed the effect of the enhanced luminosity of 
cool cores on the number of CC objects using numerical 
simulations. The authors estimate that a true fraction of 
31% strong CC objects is sufficient to explain the value of 
44% measured in the HIFLUGCS sample. Comparing the 
HIFLUGC S clusters with a sampl e of high-redshift objects 
(z > 0.5), iVikhlinin et 



al. 



(2007) noted a decrease of the 
CC fraction with redshift, which they interpret as evidence 
for a higher rate of merging events. However, thi s result 



might have been affected by a bias against CCs (Santos 
etaLl[20T0 ) 



In the following, we present our analysis of the bias in- 
troduced in X-ray cluster samples by the sharp surface- 
brightness profiles of CC clusters (hereafter, CC bias), ex- 
panding on the work of H10, with the aim of providing a 
reliable measurement of the fraction of CC vs NCC clusters. 
The reader should note that the effect investigated here and 
the one analyzed by H10 are slightly different. While H10 



analyzed the intrinsic effect induced by the increased lu- 
minosity of CC clusters, in this paper we investigate the 
bias introduced by the different detection efficiency of clus- 
ters for different surface-brightness profiles. For this work, 
we use both numerical simulations (Sect. 2|3) and obser- 
As a result, we compute the fraction 
ibjects in the local Universe taking into 
In a companion paper, we use 



vations (Sect 
of CC vs NC(To" 
account the selection bias 
our simulation tool to compute the bias expected in several 
other well-known samples and make p redictions for the fu- 
ture X-ray s urvey missions, eROS ITA ( |Predehl et al.|2010[ ) 
and WFXT ( |Giacconi et "aLl|2009| . 



Throughout the paper, we assume a ACDM cosmology 
with TO = 0.27, fl A = 0.73 and H = 72 km s _1 Mpc -1 . 



2. The simulation 

To estimate the CC bias introduced into a sample, we use a 
Monte Carlo approach where clusters are distributed follow- 
ing the observed luminosity function, with a known input 
fraction of CC vs NCC objects, and then selected using the 
criteria of the corresponding sample. A realistic description 
of the survey properties (RASS) and of the instrumental 
characteristics are implemented to perform the selection. 
The fraction of CC vs NCC clusters is then computed from 
the selected sample and compared with the input ratio. The 
dependence of the bias on several input parameters (lumi- 
nosity, redshift, absorption) is also analyzed. 

2.1. Input distributions of luminosity, redshift and column 
density 

The first step for the simulation is to select randomly a 
luminosity in the 0.1-2.4 keV band (Lx) and a redshift (z), 
following the observed distribution of clusters, i.e. 



N cl (L x ,z)=F(L 



x 



z) (-Tp) dIj x dz, 



(1) 



where the X-ray Luminosity Function (XLF) F(Lx, z) has 
the form 



F(L x ,z) = C(z)L x a eM-Lx/Mz)). 



(2) 



The dependence of the luminosity function on the redshift 
is given by L*(z) = L*(0)(l+z) A and C(z) = C(0)(l + z) B , 
where C (0 ) is the normalizatio n of the luminosity function 
at z — (Mullis et al. 20041). For the parameters of the 



luminosity function we use the values extracted from the 
BCS survey: a = 1.8 5, L*(0) = 2.8 x 10 44 hj 2 ergs s 



(Ebeling et al. 1997). The cosmological e volution of the 
:CLF follows the work of |Mullis etaL| ( |2004| ), A = -2.3 and 
B = 1.3. 

In the case of HIFLUGCS, we restrict ourselves to 
L x > 2 x 10 42 ergs s _1 and z < 0.25 to cover all the 
parameter space of the sample. At low redshift, the cosmo- 
logical dependence of the XLF is negligible, and therefore 
the redshift distribution goes like z 2 , given that the number 
of objects is proportional to <~ z 2 . The luminosity dis- 
tribution follows approximately a cut-off power law as ex- 
pected from Eq. [2j From the luminosity, we use the Lx — T 
relation, 



Lx = A 



hT ■ 
6 keV 



(3) 



Eckert, D. et al.: The Cool-Core Bias in X-ray Galaxy Cluster Samples I 



3 



to get the corresponding temperature. For the parameters 
of the relat ion, we use the values derived by |Andersson| 



et al. 



(2009), A = 1.88 x 10 45 ergs s^ 1 (2-10 keV band) and 



79. We neglect the scat ter of the relation. For gal axy 
groups, the relation is steeper ( Helsdon fe Ponman|2000[ ), so 
at low temperatures (kT < 1 keV) we use the corresponding 
determination of the Lx — T relation (A — 9.55 x 10 42 ergs 
s _1 , a = 4.9). For each value of Lx and z, we derive the 
0.1-2.4 keV flux of the corresponding object in the observer 
frame. 

In addition to the luminosity and redshift, we also sim- 
ulated a distribution of galactic hydrogen column density 
(Nh). For this purpo se, we used the 21 cm maps from 



Kalberla et al. ( 2005 1 and computed the distribution of 
Nh, restricting to galactic latitudes \l\ > 20° to match the 
considered samples. Values of Nh were then simulated fol- 
lowing the observed distribution. 



2.2. Input surface brightness profiles 

For the surface-brightness profiles of the simulated clus- 
ters, we use two different sets of parameters describing dif- 
ferent physical cases. For simplicity, we divide the popu- 
lation of clusters strictly into CC and NCC clusters (i.e. 
we neglect intermediate objects). Namely, for NCC clus- 
ters we assu me a surface-brightness profile de scribed by a 
beta model ( Cavaliere fc Fusco-Femiano||1976 1, 

(4) 



S(r) = S Q [1 + (r/r c f 



I -3/3+1/2 



while for CC clusters we assume a double-beta model, 

-3/9+1/2 



S(r) = Si ([l + (r/r cl ) 2 
+R [1 + (r/r c2 ) 2 



-3/3+1/2 



(5) 



where we use the same value of j3 for the 2 beta components, 
and R is the ratio between the two components at r = 0. 

In the first case, we choose fixed values for the rele- 
vant parameters, which were estimate d by averag i ng th e 



observed quantities from the sample of Mohr et al. ( 1999 1 
In both cases, we use a common 8 value of 0.64. For NCC 
clusters, a core radius r c — 230 kpc was chosen. For CC 
clusters, we use a core radius for the smaller beta compo- 
nent, r C 2, of 40 kpc, while the core radius of the broader 
beta component, r cl , is set to 170 kpc, a smaller value than 
for NCC clusters. All these parameters show relatively small 



scatter in the Mohr et al. (|1999|) sample. The most uncer 
tain parameter of the sinm 



bet a components at r = 
the iMohr et al.l (1 19991 



ation is the ratio between the 2 
for CC clusters, namely R. From 
sample, this quantity has a mean 
value of ~ 15 with large scatter, the values ranging from 
1.5 up to ~100 from an object to another. In the following, 
we use a constant value R = 15 for CC clusters. Indeed, 
we found that the results of the simulation do not differ 
significantly when a single value or a distribution of values 
is used for the relevant parameters (see Sect. [3]). For the 
remainder of the paper, we refer to this set of parameters 
as case I. 

As an alternative approach, we also performed simula- 
tions in which the various core radii are defined as a fixed 
fraction of r50(Q hereafter case II. More specifically, for each 

1 For a given over-density factor A, ta is defined as the radius 
within which the total mean density is a factor A above the 
critical density. 



simula ted cluster we use the scaling relation of |Arnaud et al.| 
(2005) to compute the expected value of rsoo, and choose 
for the various core radii a fixed fraction of rsoo. I n terms 
of r 50 o, for NCC systems we use r c = 0.19r 50 o, while for 
CC objects we choose r c i = 0.13r5oo and r C 2 = 0.03r5oo- 

It is unclear which of the two parameter sets should 
provide a better description of the data. Indeed, the first 
approach assumes that in cluster cores feedback effects are 
dominant, thus creating universal core properties which do 
not depend on cluster mass. Conversely, the second ap- 
proach neglects the feedback effects and assumes that the 
core properties are only determined by gravitational pro- 
cesses, in which case they should follow a self-similar rela- 
tion. Since it is known that both effects are important, our 
two choices of scale radii should provide limiting cases. 

2.3. ROSAT specifics and source detection 

The ROSAT all-sky survey (RASS) covered the whole sky 
with a typical exposu re time of a few hundred seconds 
(Cruddace et al. 2002). The instrument featured an effec- 
tive area A ~ 400 cm z @ 1 keV and a spatial resolution of ~ 
20 arcsec FWHM averaged over the FOV. The background 
count rate, which we estimated by computing the mean 
number of counts in RASS images, is low, ~ 10~ 3 counts 
s _1 arcmin~ 2 (0.1-2.4 keV). To convert fl uxes to RASS 



count rates, we used the XSPEC package (Arnaud 1996) 
vl2.6.0 and the ROSAT /PSPC response. We assumedtnat 
the spectrum of each source is describ ed by an absorbed 
MEKAL model ( [Kaastra fc Mewe||2000[ ) with varying tem- 
perature and hydrogen column density. The conversion fac- 
tor from flux to count rate was then computed for a range of 
values of kT and Nh- For any initial value of temperature 
and Nh , we estimate the conversion factor by interpolating 
the nearest points. 

In the case of extended sources, the choice of the inte- 
gration radius from the surface-brightness peak is crucial. 
Beyond a given radius, the surface brightness of the source 
becomes smaller than that of the background, and therefore 
the observed flux of a cluster in a real observation is always 
smaller than the total flux of the cluster. For this reason, 
the fluxes computed as described above cannot reproduce 
accurately the RASS. 

To overcome this difficulty, for each cluster with a given 



S(r) 2ttt dr. 



(6) 



Obviously, the surface-brightness profile S(r) is normalized 
in such a way that 



S(r) 2ttt dr = F. 



(7) 



For a cluster to match the selection criterion of HIFLUGCS 
we impose that F^g, not F, be above the flux limit of the 
sample (i.e. 2.0 x 10~ n ergs cm -2 s _1 ). 

Clearly, a realistic definition of r max is important to re- 
produce accurately the observed samples. Given that the 
original fluxes were computed insid e the maximum radius 
up to which the source is detected (Reiprich & Bohringer 



20021), we chose to use for r max the radius for which the 
surface-brightness of the source is equal to that of the back- 
ground. This should be regarded as a conservative choice, 
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Table 1. Results of different simulations of 10 6 objects with 
an input CC fraction Ci np = 0.38. Left: Surfac e-br ightness 
parameters used for the simulation (see Sect. 2.2). Right: 
Observed fraction of CC clusters in the simulated samples. 



Simulation 



Cobs 



All fixed 0.49 ± 0.02 

Distribution of /3 0.47 ± 0.02 

Distribution of r c , r c \ and r C 2 0.50 ± 0.02 

Distribution of R 0.51 ± 0.02 



10 Z 

Radius [arcsec] 



Fig. 1. Differential number-of-counts profiles (in units of 
counts per arcmin) for a CC (black) and a NCC cluster 
(blue) with the same redshift (z = 0.05) and luminosity 
(L x = 10 44 ergs s" 1 ), compared to the ROSAT /PSPC 
background profile (red). 





1<T 10" 

L x [ergs s' 1 ] 



10" 

Redshift 



10" 1 



Fig. 2. Distribution of luminosities (left) and redshifts 
(right) of the objects in the simulated flux-limited sample. 



since in most cases the source contribution can be mea- 
sured up to a somewhat larger radius. In Fig. [T] we show 
the differential number-of-counts profiles, i.e. the surface 
brightness times 27rr, for the ROSAT /PSPC background 
and two clusters, a CC and a NCC, both at z = 0.05 and 
with a luminosity of 10 44 ergs s _1 (total flux F ~ 2 x 10~ n 

ergs cnr 2 s _1 ). The value of r max in the two cases is given 3. Simulation results 



Once the simulated sample is selected, wc infer the frac- 
tion of CC clusters among the selected objects (c ovs ), and 
compare this number with the input fraction Ci np . 



by the intersection between the black/blue curve and the 



red curve. In the specific case shown here, 



490" for 



the CC case and r max = 640" for the NCC case. In terms 
of flux, this results in a ~ 25% higher observed flux for the 
CC cluster compared to the NCC cluster, for the same total 
flux. 

2.4. Construction of a simulated sample 

The critical parameter for the simulation, c, is the fraction 
of clusters hosting a cool core. Since we expect that the ob- 
served value of c in the samples is biased, we wish to com- 
pute the input value of c for which the fraction of CC clus- 
ters in a simulated sample reproduces the observed value 
in a given observed sample. Therefore, for each simulated 
cluster we randomly choose whether it is CC or NCC, im- 
posing a fraction a np of CC clusters. We then compute the 
corresponding integration radius r max as explained above 
and calculate the observed flux using Eq. [6] For a cluster to 
match the selection criterion of a sample, we then impose 
the corresponding flux limit. In addition, we also check the 
detection level of each source. More specifically, we impose 
that the signal-to-noise of a source, defined as 



We started by simulating a population according to case 
I (see Sect. 2.2). In a simulation of 10 6 objects, we found 
that ~ 900 match the selection criteria, which allows us to 
compute the observed fraction c b s with good accuracy. We 
then performed simulations for a range of values of c and 



S/N = 



N„, 



N, 



bky 



(8) 



where N source and N^kg are the total number of counts 
of the source and the background integrated up to r max , 
be above 5.0 to be selected in the sample. In the case of 
HIFLUGCS, all the objects are detected with S/N > 30, 
so this condition is unnecessary, but in general this is not 
true. 



found that the observed value of 49% ( Chen et al.||2007 ) 
corresponds to an input ratio Ci np = 0.38 ± 0.02, i.e. to a 
significant bias of ^29%. In case II (see Sect. 2.2), we find 
a similar result (ci np — 0.37 ± 0.02). We also investigated 
how the bias changes when we use distributions of the rel- 
evant parameters rather than fixed values. The results are 
shown in Table [TJ using the same input value ci np = 0.38. 
Introducing distributions of parameters has little influence 
on the output results, so for the remaining of the paper we 
use fixed values for the surface-brightness parameters. 

In addition to the total bias, it is also interesting to in- 
vestigate the dependence of c ous on other quantities, in par- 
ticular on luminosity, redshift and hydrogen column den- 
sity. A simulation with a very large number of objects is 
required to study the dependence of c {, s on these param- 
eters. Fixing the value of Ci np to 0.38, we then performed 
large simulations of 10 s objects for the two parameter sets 
(case I and case II), where as expected ~ 90,000 match 
the selection criterion. Figure 2] shows the luminosity and 
redshift distributions of the selected objects. The distribu- 
tions peak at a luminosity ~ 3 x 10 44 ergs s _1 and a red- 
shift ~ 0.05, which matche s well the observed quantitie s 
of the HIFLUGCS sample ( |Reiprich fc Bohringer| |2002l . 
Moreover, as expected, the logAMogA' of the sample is well- 
represented by a power law with an index —3/2, in agree- 
ment with the properties of the actual sample. 

We studied the dependence of c oos on redshift, luminos- 
ity and absorption. The resulting plots are shown in Fig. [3] 
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Fig. 3. Dependence of the observed fraction of CC objects, c fe s , on the 0.1-2.4 keV luminosity (left) and the redshift 
(right). In both cases, the black or blue data show the simulation results for case I or case II, respectively (see text). The 
dashed red lines show the input fraction of CC objects, Ci np , which was fixed to 0.38 for the simulations. 



(black: case I; blue: case II). Since they represent extreme 
cases where the core radii are determined only by gravi- 
tation (blue) and feedback (black) effects, the two curves 
should give lower and upper bounds to the actual values. In 
case I, we observe a strong dependence of c f, s on the lumi- 
nosity. In particular, the low-luminosity objects (Lx < 10 43 
ergs s _1 ) are very strongly biased towards CC clusters. 
Since the number of objects in this luminosity range is rela- 
tively small (see the left panel of Fig. [2]), this effect does not 
have a big influence on the mean value of c b s - This effect 
is related to the background level. Indeed, these objects are 
both low-luminosity and nearby, so their surface-brightness 
peak is small. In the case of some NCC clusters, it can 
even be below the surface brightness of the background, in 
which case r max = and the source is not detected at all. 
Conversely, because of their peaked profiles CC objects with 
the same total luminosity are easily detected. Therefore, 
CC clusters completely dominate the detected objects at 
luminosities below 10 43 ergs s _1 . In case II, this effect, al- 
beit present, is much less severe, because low-luminosity 
(i.e. low-mass) systems have smaller core radii. As a re- 
sult, they are more concentrated, and t he effect becomes 
less important. It has been claimed (e.g., Chen et aL||2007 ) 
that low-mass systems (groups, poor clusters) are predom- 
inantly CC. If the CC bias is not taken into account, we 
see from our simulation that one would immediately reach 
such a conclusion, even if the true fraction of CC objects 
would not depend at all on luminosity. 

From Fig. [3] we also observe a dependence of c b s on the 
redshift of the object. In particular, we note a significant in- 
crease in c ft s . at the highest rcdshifts. This is a consequence 
of the cut-off in the luminosity function. At rather high 
redshifts, only the most luminous objects are selected, and 
above the cut-off luminosity their number decreases quickly 
with increasing luminosity. Since, for a given limiting flux, 
the associated luminosity is smaller for CC than for NCC 
systems, the relative number of CC objects will be signifi- 
cantly larger. We note that unlike the low-luminosity effect, 
which is due to a combination of the intrinsic properties of 
the different population and of the background level, this 
effect is caused by the flux limit. With a lower flux limit, 
this effect would not disappear, but would be shifted to 



higher redshift, where the effect of the cut-off in the lumi- 
nosity function will start to appear. Unlike the luminosity 
dependence, in case II we find a stronger evolution of c b s - 
In this case, since we are only selecting high-mass systems 
the core radii are proportionally larger, which increases the 
bias. 

As expected, the selected sample shows very little de- 
pendence on Njj- Indeed, the measured value of c f, s is con- 
stant for varying Njj, and the total fraction of CC is very 
similar to the one obtained when the effect of absorption is 
completely neglected. 

To visualize the simulated objects and compare them 
with the actual data, we plotted all the detected clusters 
in the Lx — z plane, and ad ded to the p l ot the data of 
the HIFLUGCS sample from IChen et all ((20071). In this 



case, we classified the clusters as CC/WCC on the basis of 
their classical mass deposition rate, CC clusters being the 
ones showing a non-zero mass deposition rate. The resulting 
plot for case I can be seen in Fig. |2J Black (NCC) and red 
(CC) dots each represent a single simulated cluster, while 
the green symbols (NCC) and yellow circles (CC) show the 
actual data for the objects of the extended HIFLUGCS 
sample. 

Interestingly, we see in Fig. [3] that at all redshifts the 
simulated CC clusters populate a broader range of lumi- 
nosities than NCC clusters. In other terms, the total flux 
limit of the sample is different for CC and NCC clusters. 
At low luminosities (\ogLx < 43.5) CC objects dominate 
the observed and simulated data, which again indicates an 
observational bias (see the left panel of Fig. [3]) . All but one 
NCC clusters lie in the area populated by the black dots, 
while the remaining object is found just below the limit. 
This indicates that our simulation reproduces well the prop- 
erties of the sample. A similar conclusion is reached for case 
II. 



4. X MM- Newton and ROSAT analysis of 

surface-brightness profiles from the HIFLUGCS 
sample 

We analyzed available data for the HIFLUGCS sample, 
with the aim of providing an alternative estimate of the 
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Fig. 4. Comparison between simulated objects and observational data in the log Lx vs z plane. The black and red dots 
represent the simulated NCC and CC clusters, respectively. The green symbols (NCC) and yellow circles (CC) show 
the data from the extended HIFLUGCS sample ( Chen et al.||2007 |. The dashed blue line shows the flux limit for point 
sources. 



CC bias and of extracting a subsample that is unaffected by 
the CC bias. We proceeded with the systematic analysis of 
surface-brightness profiles for all clusters in the HIFLUGCS 
sample. Since the sample is mostly composed of nearby 
objects, the Field Of View (FOV) of the selected instru- 
ments is a key feature to observe regions as large as pos- 
sible. Indeed, we are interested in the extraction of total 
fluxes within a given physical radius, so a detailed charac- 
terization of the central regions is not necessary, but it is 
important to detect the emission to relatively large radii. 
As a result, we decided to use XMM-Newton /EPIC and 
ROSAT /PSPC pointed data to pursue our goals, since they 
are the X-ray instruments with the largest FOV. In most 
cases, when both XMM-Newton and ROSAT pointed data 
were available, we combined the surface-brightness profiles 
from the two instruments to take advantage of the better 
XMM-Newton point spread function (6 arcsec FWHM) in 
the central regions and of the larger FOV (1 square degree) 
and lower background of ROSAT in the external regions. 
For very nearby objects (z < 0.015), the ROSAT PSF is 
sufficient to resolve the cores, and a detection to large radii 
is needed to constrain the parameters, so we restricted our 
analysis to ROSAT. Conversely, for the more distant ob- 
jects (z > 0.08) the apparent size of the objects is smaller, 



so a better PSF is necessary to resolve the core. In these 
cases, only XMM-Newton was used. 



To construct a subsample free of selection bias, we mea- 
sure fluxes in an annulus around the core to exclude the 
core component. We then perform a new selection on the 
basis of these fluxes. Generally, the flux of a cluster is not a 
well-defined quantity, since it is computed as the total flux 
within the area where the source is detected. However, in 
this work we wish to extract fluxes within a well-defined 
physical region, which we choose to be t^qq. Since in most 
cases we do not detect the sources at this radius, a best-fit 
surface-brightness model is used to e xtrapolate the fluxes 
to r 50 o- We use the scaling relations of |Arnaud et al. (2005) 
and cluster virial temperatures from H10 to estimate r$oo 
for each object. Overall, this corresponds to a size which 
ranges from ~ 500 kpc for small, group-like objects up to 
~ 1,500 kpc for the most massive, hot clusters. The val- 
ues (physical and apparent) of r 50 o, the redshift and virial 
temperatures of all HIFLUGCS objects are summarized in 
Table |1J together with the log of the available data. 
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4.1. Data analysis 
4.1.1. XMM-Newton 

XMM-Newton /EPIC data analysis was carried out in a sys- 
tematic way using the standard XMMSAS software v9.0. 
For each EPIC detector, light curves in the soft (2-5 keV) 
and hard band (10-12 keV) were extracted, and events were 
filtered to exclude flares. To minimize the instrumental 
background, we extracted images fro m the cleaned even t 



lists in a narrow band (0.7-1.2 keV, Ettori et al. 
which excludes the prominent Al 



2010) 
1.5 keV) ancT gT 



1.8 keV) background lines (Leccardi & Molcndi 2008J). As a 
result, the total background is dominated by the sky back- 
ground, which is convolved with the vignetting of the in- 
strument in the same way as the source signal. For this 
work, since we are interested in the characterization of the 
total flux it is important to detect the sources to large radii. 
Given that the objects considered here are bright, combin- 
ing the three EPIC detectors is not necessary. We chose to 
use MOS because of its better imaging capabilities, and for 
simplicity we restricted our analysis to MOS2. To extract 
surface-brightness profiles, we generated MOS2 images with 
4 arcsec spatial binning in the 0.7-1.2 keV band, and the 
corresponding exposure maps using the SAS tool eexpmap. 

4.1.2. ROSAT 

When available, we restricted our analysis to 
ROSAT/PSPC pointed observations. We used the fil- 
tered event files available in the HEASARC archive^ and 
extracted images using the XSELECT tool available within 
the HEASOFlj^] software package with 15 arcsec bins in 
the 0.4-2.0 keV bands. Exposure maps in the same energy 
band were generated using the PSPC detector maps and 
the instrument attitude files through the pcexpmap tool. 

4.2. Construction of surface-brightness profiles and flux 
measurements 

For every cluster and instrument, we used the output 
count images and corresponding exposure maps to extract 
surface-brightness profiles. After the detection of individual 
point sources using a local background method, we excised 
from the images the corresponding regions. We then ex- 
tracted total count profiles, centered on the image centroid, 
defined as 



«*),<</» = 



X/ Ix,y X/ I%,y 



(9) 



where I XyV is the count rate in the pixel with image coordi- 
nates (x,y). The minimum bin size of the profile was set to 
7 arcsec (XMM-Newton) and 20 arcsec (ROSAT), and the 
bins were grouped to ensure a minimum of 50 counts per 
bin and permit the use of x 2 statistics. The count profiles 
were then divided by the corresponding exposure to correct 
for the instrument vignetting. 

As a result, we obtained exposure-corrected profiles, 
which we fitted by a source 4- background model. A double- 
beta profile was preferred to the standard beta profile when 
the improvement to the fit was found to be statistically sig- 
nificant by the F-test. The background-subtracted source 



profile SB(r) was then inferred. To ensure that the source 
dominates over the background, this procedure was applied 
up to the radius r max where the model source intensity is at 
least 2 times that of the background. The best-fit surface- 
brightness model was then used to extrapolate the fluxes 
to rsoo- We define the total fluxes within rsoo as: 



F tot = SB(r)w(r) 



S(r) 2nr dr, 



(10) 



2 http://heasarc.gsfc.nasa.gov/ 

3 http: / /heasarc. nasa.gov/lheasoft / 



where SB(r) represents the measured background- 
subtracted surface-brightness profile in the bin correspond- 
ing to the radius r, w(r) is the width of the bin, and S(r) 
is the best-fit surface-brightness model. In the radial range 
where we use the actual data, the error on the measured 
flux is simply calculated by propagating the Poisson statis- 
tics. In the outer regions where we use the model surface- 
brightness profile to estimate the flux up to rsoo, we esti- 
mate the error through a Monte Carlo approach: randomly- 
picking the model parameters, we compute the correspond- 
ing flux for each parameter set, and extract the resulting 
flux distribution. Since the core contribution is expected to 
be negligible in this range, when a double-beta model is re- 
quired we fix the parameters of the inner beta component. 
Moreover, the outer core radius r c \ and the /3 parameters 
are strongly degenerate, so we only let /3 vary within its un- 
certainties. The error on the second term of Eq. [10] is then 
given by the RMS of the resulting flux distribution. Overall, 
in the vast majority of cases the statistical uncertainties on 
the extracted fluxes are small (F/AF ~ 50), however in 
some cases where the parameters are poorly constrained 
the relative uncertainty on the flux measurement from to 
r5oo can be quite large. 

In a majority of cases (38 objects out of 64), we com- 
bined ROSAT and XMM-Newton surface-brightness pro- 
files and fitted them jointly to get better constraints on 
the parameters. As an example, Fig. [5] shows the XMM- 
Newton /MOS2 (blue) and ROSAT/PSPC (red) surface- 
brightness profiles of the cluster A3571 fitted by a double- 
beta profile plus a constant for the background. While the 
physical parameters of the model where fitted jointly, the 
individual normalizations and background levels where ad- 
justed independently. In the innermost 1 arcmin of the pro- 
file, ROSAT data were ignored to avoid problems of mixing 
caused by the broader PSF. 

To convert the count rates into physical fluxes, we used 
the XMM-Newton /MOS2 and ROSAT/PSPC redistribu- 
tion matrices and on-axis effective area, and used XSPEC 
to fold the spectrum of each source, modeled by an ab- 
sorbed MEKAL model with the appropriate temperature 
and Njj, with the appropriate instrument response. The 
folded spectra were then used to convert the extracted 
count rates into unabsorbed physical fluxes in the 0.5-2.0 
keV band. 

As a byproduct of this work, we were able to esti- 
mate the cross-calibration between the two instruments 
(see Appendix |A| for details). Indeed, since they are bright, 
persistent sources, clusters of galaxies are ideal sources to 
cross-calibrate between two independent missions. On aver- 
age, we find that MOS2 gives ^15% higher fluxes compared 
to PSPC, and hence there is a clear cross-calibration issue 
which must be taken into account. Since the HIFLUGCS 
sample was selected using fluxes extracted from ROSAT 
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Fig. 5. XMM-Newton/MOS2 (blue) and ROSAT /PSPC 
(red) surface-brightness profile of the cluster A3571, fitted 
jointly with a double-beta model. The individual normal- 
izations and background levels were fitted independently. 
The higher XMM-Newton data points and model at large 
radii are caused by a higher background. The bottom panel 
shows the deviations (in a) from the surface-brightness 
model. 



only, we decided to use the ROSAT normalization to ex- 
tract fluxes, and apply a 15% cross-calibration factor to 
the XMM-Newton fluxes. After applying such a renormal- 
ization, XMM-Newton and ROSAT fluxes agree with -5% 
scatter. 

The best-fit parameters and 1-a errors, as well as the 
total fluxes from XMM-Newton and ROSAT are summa- 
rized in Table [5} The flux limit of the HIFLUGCS sample 
(2.0 x 10 -11 ergs cm -2 s _1 in the 0.1-2.4 keV band) trans- 
lates into a total flux of 1.2 x 10~ n ergs cm -2 s _1 in the 
0.5-2.0 keV band for a typical temperature of 4 keV. All 
XMM-Newton fluxes should be rescaled by 15% when com- 
pared to ROSAT to account for the systematic calibration 
difference. Figure [6] shows the distribution of ratio between 
the ROSAT fluxes ex tracted following Eq. 10 and t he orig- 
inal fluxes quoted by |Reiprich fc Bohring er (2002). In the 
vast majority of cases, we can see that the fluxes estimated 
using Eq. |10| agree very well with the original fluxes. In three 
cases, our fluxes differ from the original fluxes by a fac- 
tor ~ 2, either higher (Fornax and A1367) or lower (NGC 
4636). All of these are nearby, very extended objects. The 
difference can be explained by the low surface brightness of 
the outer regions, which prevents the detection of a large 
part of the flux, and by the instability of the extrapola- 
tion over a large radial range. Overall, this is however not 
a problem, given that our definition of cluster fluxes differs 
from the original method, where the fluxes were measured 
by integrating the flux within the area where the clusters 
are detected. 

According to our analysis, one cluster (Zw III 54) does 
not have a sufficient total flux to match the selection crite- 
rion of the sample. This cluster is the weakest of the sample 
and was at the extreme limit of the original selection, and 
therefore, even though our redetermined flux does not dif- 
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Fig. 6. Distribution of the ratio between our redetermined 
fluxes and the original fluxes from Rciprich & Bohringer 
(2002). Since the original fluxes were computed in the 0.1- 
2.4 keV band, they were rescaled to match our choice of 
energy band (0.5-2.0 keV). The most prominent outliers 
are highlighted. 



fer significantly from the original one, we find its flux to be 
below the flux cut. 



5. The unbiased subsample 

5.1. Selection of the subsample 

To measure the CC fraction, we wish to extract a sub- 
sample of HIFLUGCS which should be free of the CC 
bias. The fraction of CC objects would then be computed 
from the subsample and not from HIFLUGCS. At variance 
with the cores, it has been found that in cluster outskirts 
the surface-brightness profiles are essentially self-simila r 
( |Neumann||2005[ iCroston et aT|2008| |Leccardi et~aL|2010[ ) , 



so the selection based on the fluxes in the outer regions 
should be unbiased. We define a new sample by selecting 
those objects for which the flux in an annulus excluding 
the core, and not the total flux, is above a given flux limit. 
While the choice of the outer radius is straightforward (we 
select rsoo), the same is not true for the inner radius, Ti n . 
Here we require a value such that the core emission would 
be excluded as much as possible, but that would still allow 
a statistically accurate determination of the fluxes. To this 
end, we analyzed the distributions of the best-fit parame- 
ters from Table O In the cases where a double-beta model 
was required by the fit, we searched for a typical radius 
at which the contribution of the core can be considered to 
be negligible. To do this, we looked for the radius above 
which the flux can be attributed mostly to the outer beta 
component. All in all, we find that ri n = 0.2r5oo is a reason- 
able choice. With such a choice, in average 97% of the flux 
is attributed to the main beta component, and the fluxes 
can still be computed with good accuracy. For comparison, 
this radius is larger than or comparable to the cooling ra- 
dius r coo i f or whi ch the cooling time exceeds 7.7 Gyr (H10). 
Maughan (2007) excised the region within 0.15r 50 o to ex- 
clude the effects of the core. In this case, in average ~ 6% 
of core flux is still included. For the remaining of the paper, 
we compute the fluxes in the 0.2rsoo — ?~5oo radial range. 
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To select the corresponding flux limit in this radial 
range, we considered a beta profile with the mean parame- 
ters from our analysis and computed the fraction of the flux 
integrated in this range. From Table [5j we find /3 ~ 0.64 
with a scatter of 0.13 and r c \ ~ 0.16rsoo with a scatter of 
0.08r 50 o- Integrating a surface-brightness profile with such 
parameters, we estimate that 59% of the flux originates 
from the 0.2rsoo — f^oo, an d therefore the HIFLUGCS flux 
limit translates into a minimum flux of 7.2 x 10~ 12 ergs 
cm~ 2 s" 1 (0.5-2.0 keV band, 0.2rsoo — rsoo)- Because of 
the 15% normalization difference, this corresponds to a flux 
limit of 8.3 x 10~ 12 ergs cm~ 2 s _1 for XMM-Newton fluxes. 



5.2, Incompleteness 

Since our sample is selected as a subset of HIFLUGCS and 
not directly from RASS data, there may be objects that 
are not in the HIFLUGCS sample but should be in ours. 
In other terms, there could be some objects with a flat 
profile, for which the total flux would be slightly below the 
HIFLUGCS limit but the flux in the 0.2r5oo — ?"5oo radial 
range would be above our new limit. This effect could have 
an influence on our measurement of the CC fraction using 
our subsample, and it must be quantified. 

To estimate the number of objects which might be 
missed in our subsample, we computed the distribution of 
the ratio between our fluxes in the 0.2rsoo — ?"5oo and the 
original survey fluxes, excluding the most prominent out- 
liers from Fig. [6] and rescaling the distribution by the ratio 
of the two flux limits. This defines a probability distribution 
PDF(x) for the probability that the flux of a cluster would 
be increased or decreased with respect to the flux limit 
when going from the original selection in the HIFLUGCS 
sample to the one we apply. In the probability distribution, 
the peaked objects are the ones with a ratio much below 1, 
while the flat systems have a ratio larger than 1. In Fig. [7] 
we show the resulting probability distribution. 

For an object with a given total flux F < F[ im , the 



probability to have a 0.2rsoo 
flux limit is given by 



/>oo 

P(F) = / PDF 

JF Hm /F 



(x) dx, 



rsoo flux greater than our 



(11) 



We can then estimate the number of missed objects by inte- 
grating the probability distribution weighted by the number 
of objects at a given flux, i.e. 



^missed — ^ ^ 



N obj (F) x P(F). 



(12) 



where N b j is estimated from the logiV — lo g S of the 
REFLEX ([Bohringer et al.|2004| and NORASS |B5hringer 
et al. 2000) samples from which HIFLUGCS is extracted. 
As a result, the average number of objects which should 
be selected in our subsample but not in HIFLUGCS is 
Nmissed ~ 3.89. Given that the objects which might be 
missed are the ones with a high flux outside the core, their 
surface-brightness profile should be very flat, so we expect 
them to be mostly NCC. In conclusion, we expect that sta- 
tistically between 3 and 4 NCC clusters are missed from 
our selection. 
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Fig.7 

f' 

computed in the 0.2r5oo 

responding flux limit (7.2 x 10~ 12 e rgs cm~ 2 s" 1 , 0.5-2.0 
keV), and f id is the original flux from |Reiprich k, Bohring er| 
( 2002 ) rescaled by the corresponding flux limit (2.0 x 10" 
ergs cm -2 s -1 , 0.1-2.4 keV). This defines a probability dis- 
tribution PDF{x) for the probability that the flux of a 
cluster is increased or decreased when performing our new 
selection. 



5.3. Estimate of the CC fraction 

After the selection of the clusters based on their flux in 
the 0.2r5oo — f*5oo radial range, we reject 13 clusters from 
the original sample. Table [2] shows the corresponding ob- 
jects and fluxes, along with their classification by H10. 
Interestingly, all the rejected objects were classified as CC: 9 
of the clusters were classified as strong CC clusters (i.e. clus- 
ters with a central cooling time < 1 Gyr, hereafter SCC), 
and the remaining 4 as weak CC (central cooling time be- 
tween 1 and 7.7 Gyr, hereafter WCC). Among the rejected 
WCC clusters, two of them (A2244 and A1650) have a cen- 
tral cooling time just above 1 Gyr, so they are borderline 
objects which resemble closely strong CC objects. One ob- 
ject, Zw III 54, has a redetermined total flux below the 
HIFLUGCS cut (see above). As a result, our selection only 
affected objects which exhibit strong CC characteristics and 
a prominent central emission excess. We also note that an 
important fraction of the clusters (5 out of 13, namely EXO 
0422, Zw III 54, A 3581, NGC 4636, and Sersic 159-03) have 
a temperature below 3 keV, i.e. they are low-mass, low- 
luminosity objects. This is in qualitative agreement with 
our simulation results (see Fig. |3|, which show a strong 
bias of the original sample in the case of low-luminosity ob- 
jects. Moreover, 2 of the remaining objects (RX J1504 and 
A2204) are distant (z > 0.15), strong CC objects, again in 
agreement with the predictions of our simulations (see the 
right panel of Fig. [3]) . 

After the selection, our subsample comprises 51 objects 
(see Table [3]). The table shows our selection of clusters, 
together with their fluxes in the 0.2rsoo — ?"5oo annulus, 
central cooling time and classification from the work of H10. 
Four clusters (MKW 4, MKW 3S, A2163 and A2589) show 
a flux below the limit, but consistent with it within lcr, so 
we included them in our selection. 
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Table 2. HIFLUGCS clusters rejected for the present work. 
Column description: 1: Cluster name; 2 and 3: 0.5-2.0 
keV fluxes in the radial range 0.2rsoo — rsoo from XMM- 
Newton /MOS2 (2) and ROSAT /PSPC (3), in units of 
10 -12 ergs cm~ 2 s _1 ; 4: central cooling time in units of 



z-l/2 
"71 

(SCC 



Cyr, from H10; 5: cluster classification from H10 
= strong cool-core, WCC = weak cool-core). 



Cluster 


Fx MM 


Fro SAT 


CCT 


Class 


A133 


6.88±0.38 


5.81±0.24 


0.47 


SCC 


A3112 


7.54±0.24 


6.15±0.21 


0.37 


SCC 


Zw III 54 


4.90±0.82 




5.48 


WCC 


EXO 0422 


5.03±0.67 




0.47 


SCC 


NGC 4636 


2.78±0.24 


2.52±0.22 


0.21 


SCC 


A1650 


5.96±0.55 




1.25 


WCC 


A1651 


8.25±0.28 


6.83±0.20 


3.63 


WCC 


A3581 


5.87±0.38 




0.52 


SCC 


RX J 1504 


2.71±0.19 




0.59 


SCC 


A2204 


4.93±0.38 




0.25 


SCC 


A2244 




5.02±0.26 


1.53 


WCC 


Sersic 159-03 


4.61±0.36 


3.87±0.25 


0.88 


SCC 


A2597 


3.64±0.22 


2.88±0.18 


0.42 


SCC 



Of the 51 objects in our subsample, 19 (37%) are clas- 
sified as strong CC, 18 (35%) as NCC, and 14 (27%) are 
WCC. When taking into account the cl uste rs which might 
be missed in our subsample (see Sect. 5.2) and assuming 
that all of them are NCC, the fraction of SCC is reduced to 
35% and that of NCC increased to 40%. Compared to the 



analysis of |Mittal et aL]p009[ ), which found 44% SCC clus- 
ters, our value of 35% corresponds to a bias of 26% in the 
original sample. This result is in excellent agreement with 
the bias estimated from our simulation (29%, see Sect. [3]). 



5.4. The ratio F a 



nnuLus 



/F a 



as a cool-core indicator 



Given that all the clusters excluded by our analysis present 
CC characteristics, it is clear that the ratio between the 
flux in the core, from to 0.2r5oo (hereafter F core ) and 
the flux in the annulus 0.2r 500 - 7-500 {Fannuius) strongly 
depends on the state of the cluster. In other terms, the 
ratio F annu i us / Fcore could be used as an indicator of the 
CC state. To investigate this possibility, we checked for a 
positive correlation between this quantity and the central 
cooling time (CCT) as found by H10. In Fig. [8] we show 
the ratio F annu i us / F core as a function of the CCT for all 
64 HIFLUGCS clusters, fitted by a simple power law. A 
positive correlation is indeed observed, with the best-fit pa- 
rameters given by 

TP f CCT \ °- 557±0 - 008 

fannuius = (Q m ± Q QQg-J x ^±±_ \ (13) 



F, 



VI Cyr/ 



The correlation between the two quantities is significant: 
Spearman's correlation factor for the relation is p = 
0.73t^o5- The corresponding likelihood of the correlation 
occurring by chance is < 10~ 5 . The scatter of the rela- 
tion is ~ 24%. Given that the CCT is one of the most 
commonly-used CC indicators, we conclude that the quan- 
tity F annu i us / F core can indeed be used as an indicator of 
CC type, and thus Eq. [13] can be used to estimate the cen- 
tral cooling time. This indicator is robust, since it does 
not depend on any d eprojected quant it y (sim ilar to the 
"cuspiness" indicator, Vikhlinin et al. ( 2007[ )). However, 



Table 3. Unbiased HIFLUGCS subsample used in the 
present work. Column description: 1: Cluster name; 2 and 
3: 0.5-2.0 keV fluxes in the radial range 0.2rsoo — rsoo from 
XMM-Newton/MOS2 (2) and ROSAT /PSPC (3), in units 
of 10 -12 ergs cm~ 2 s _1 ; 4: central cooling time in units 

— 1/2 

of h 71 Gyr, from H10; 5: cluster classification from H10 
(SCC = strong cool-core, WCC = weak cool-core, NCC = 
non cool-core). 



Cluster 


Fx MM 


Fro sat 


CCT 


Class 


A85 


23.45±0.73 


20.38±0.38 


0.51 


SCC 


A119 


22.17±0.78 


19.11±0.43 


14.03 


NCC 


NGC 507 


9.35±0.36 


8.32±0.34 


0.48 


SCC 


A262 


30.52±1.87 


24.62±1.16 


0.43 


SCC 


A400 


11.69±0.71 


11.36±0.31 


8.04 


NCC 


A399 


13.74±0.85 




12.13 


NCC 


A401 


17.99±0.48 


16.38±0.35 


8.81 


NCC 


Fornax 




80.88±10.79 


0.69 


SCC 


2A 0335 


16.03±1.51 


12.78±1.68 


0.31 


SCC 


A3 158 


14.18±0.39 


12.97±0.34 


8.22 


NCC 


A478 


10.57±0.20 


9.44±0.14 


0.43 


SCC 


NGC 1550 


24.27±4.96 




0.23 


SCC 


A3266 


26.80±0.41 


21.08±0.36 


7.62 


WCC 


A496 


27.42±1.37 


23.92±1.18 


0.47 


SCC 


A3376 


10.75±0.68 


10.71±0.45 


16.47 


NCC 


A3391 




8.667±0.72 


12.46 


NCC 


A3395s 


11.16±1.56 


9.955±1.23 


12.66 


NCC 


A576 


10.90±0.81 




3.62 


WCC 


A754 


31.25±0.77 


24.74±0.48 


9.53 


NCC 


Hydra A 


10.52±0.21 


8.96±0.16 


0.41 


SCC 


A1060 




25.81±4.33 


2.87 


WCC 


A1367 




40.50±1.03 


27.97 


NCC 


MKW 4 




6.12±1.86 


0.28 


SCC 


Zw CI 1215 


9.09±0.53 




10.99 


NCC 


A3526 




72.64±11.28 


0.42 


SCC 


A1644 


21.77±1.99 




0.84 


SCC 


Coma 




128.59±1.74 


15.97 


NCC 


NGC 5044 




12.57±0.48 


0.21 


SCC 


A1736 


14.65±0.75 




16.59 


NCC 


A3558 


27.38±0.26 


26.10±0.23 


1.69 


WCC 


A3562 


13.81±0.19 


11.24±0.25 


5.15 


WCC 


A3571 


38.61±2.04 


34.91±1.38 


2.13 


WCC 


A1795 


14.97±0.33 


12.14±0.22 


0.61 


SCC 


MKW 8 


9.22±4.94 




10.87 


NCC 


A2029 


16.06±2.53 


13.42±1.57 


0.53 


SCC 


A2052 


11.69±0.75 


9.55±0.61 


0.51 


SCC 


MKW 3S 


8.34±0.71 


6.67±0.53 


0.86 


SCC 


A2065 


9.63±2.03 




1.34 


WCC 


A2063 


12.38±1.57 


11.02±1.18 


2.36 


WCC 


A2142 




18.38±0.33 


1.94 


WCC 


A2147 




31.31±1.26 


17.04 


NCC 


A2163 


8.06±0.30 




9.65 


NCC 


A2199 


26.93±0.83 


22.13±0.49 


0.60 


SCC 


A2256 


26.28±0.70 


23.21±0.36 


11.56 


NCC 


A2255 


9.52±0.29 


7.86±0.14 


20.66 


NCC 


A3667 


35.86±0.35 


29.73±0.27 


6.14 


WCC 


A2589 


8.05±0.51 


6.98±0.35 


1.18 


WCC 


A2634 




11.12±0.86 


1.52 


WCC 


A2657 


9.38±0.48 


8.94±0.21 


2.68 


WCC 


A4038 


18.50±0.56 


16.31±0.45 


1.68 


WCC 


A4059 


10.00±0.38 


8.161±0.32 


0.7 


SCC 



tight constraints on the surface-brightness profile parame- 
ters are necessary, since the fluxes need to be extrapolated 
to rsoo- This is a problem for very extended, nearby objects, 
for which the uncertainties when extrapolating to rsoo are 
quite large. Indeed, all of the outliers in Fig. [8] are nearby 
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Fig. 8. Ratio between the flux in the 0.2r 50 o — ?"5oo radial 
range and the flux in the core integrated up to 0.2r5oo, 
Fannuius / F core (this work) , as a function of the central cool- 
ing time in units of Gyr (H10). The red solid line shows a 
fit to the data with a power law. Spearman's correlation 
coefficient for the relation is p = 0.73+ ; 05 . 



(z < 0.02), very extended objects, for which it is difficult to 
extrapolate the fluxes. Moreover, a number of these objects 
(NGC 507, A262, NGC 1550, MKW 4) are galaxy groups. 
Because of the difficulty of observing their outer regions, 
for groups the measured values of f3 te nd to be flatter than 
the real value ( Gastaldello et al7||2007 1 , which probably re- 
sults in an over-estimation of F annu i us . On the other hand, 
this indicator is excellent for intermediate-redshift objects 
(z ~ 0.1), for which the surface-brightness profile parame- 
ters can be well constrained, and for which PSF smearing 
can make deprojection analyses with XMM-Newton diffi- 
cult. 



6. Discussion 

Expanding on the work of H10, our analysis reveals that 
the flux- limited samples of galaxy clusters, and in partic- 
ular the HIFLUGCS complete sample, are significantly bi- 
ased in favor of CC objects, and hence this effect must be 
taken into account when estimating the fraction of CC ob- 
jects from a flux-limited sample. Our assessments of the 
bias through observations and simulations agree very well 
with each other, which strengthens our point. Overall, we 
find that ~ 29% of the objects should be removed from the 
original HIFLUGCS sample if one wishes to estimate the 
CC fraction, and we extracted a subsample as free as pos- 
sible of the bias, which can also be used for other purposes 
(see Table [3}. 

Once the appropriate corrections have been applied, we 
find that considerably less than half of the objects (38% in 
our subsample, 35% when taking incompleteness into ac- 
count; see Sect. 5.3) are classified as strong CC clusters, 



i.e. exhibit a central cooling time below 1 Gyr. This re- 
sult has important repercussions on the structure formation 
scenarios. At varia nce with the orig inal predictions of the 
cooling-flow model (Peres et al. 1998 ), in the scenario where 
CC clusters trace relaxed objects our result implies that a 
majority of clusters has not reached a stable state. In other 
terms, in the majority of objects injection of entropy in the 



ICM by merging events or giant AGN outbursts preve nts 
the formation of cool cores (McNamara & Nulsen 2007). 

Obviously, if cool cores cannot be destroyed efficiently 
by merging events and the state of a cluster is defined once 
and for all during the cluster formatio n process, as sug 



gested by some nume rical simulations (Poole et al. 2008 
McCarthy et al.l 2008 1 , the fraction of CC can only be pre 



dieted through large cosmological simulations. Conversely, 
assuming that clusters evolve through time, the CC frac- 
tion can be used to estimate roughly the rate of entropy 
injection events in the local Universe. 

Assuming that the NCC fraction evolves through time, 
the evolution of the fraction of NCC objects can be de- 
scribed as 



df. 



NCC 



— T he fee — f 'cool j NCC, 



(14) 



where r^ e is the rate of heating events per dark matter 
halo, r coo i is the cooling rate, Jncc is the fraction of both 
NCC and WCC clusters and fee is the fraction of SCC 
clusters. This equation simply reflects the idea that heat- 
ing events will transform CC clusters into NCC (or WCC), 
while cooling will do the opposite. The cooling rate is given 
by r coo i — 1/tr, where the relaxation timescale tr is the 
time needed for a cluster to relax from a major heating 
event. In this framework, SCC clusters are the ones which 
have not experienced any major entropy injection event 
during the relaxation timescale, while NCC and WCC have 
experienced or are currently experiencing an entropy injec- 
tion phase. 

If we assume that the evolution of fNCC is slow, we are 
close to a stationary situation, and Eq. [14] reads 



The 



i cool ' 



fNCC 1 fNCC 



CC 



t r fee 



(15) 



From the simulations of Poole et al. (2006), it appears that 
clusters require in average a timescale of ~ 5 Gyr to relax 
from a merging event. This also corresponds to the mean of 
the distribution of central cooling times, so we use t r = 5 
Gyr for our calculation. In the ACDM cosmology, a look- 
back time of 5 Gyr corresponds to a rcdshift of 0.4, so we 
are probing the rate of entropy injection events integrated 
up to redshift 0.4. As a result, we find that the rate of major 
heating events per cluster is roughly 



1/3 Gyr" 



(16) 



This number can be compared with the predictions 
from full cosmological simulations (see e.g., Fakhouri & 
|Ma| [2008). Obviously, this estimate relies on several strong 
assumptions, in particular on the estimate of the relax- 
ation timescale and on the approximation of equilibrium 
From a statistica l analy sis of cold fronts in galaxy clusters 



Ghizzardi et al. (2010) estimate a merger rate of ~ 1/3 
Gyr -1 , similar to ours. These numbers can however not be 
directly compared, since most of the merging events respon- 
sible for producing cold fronts are minor mergers which are 
not capable of disrupting a cool core. On the other hand, 
our calculation gives an estimate of the rate of major en- 
tropy injection events, whether they are mergers or AGN 
outbursts, which are capable of transforming a CC cluster 
into NCC or WCC. 

An interesting result of our simulation is the strong de- 
pendence of the CC bias on luminosity (see the left panel 
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of Fig. While for clusters (L x > 3 x 10 43 ergs s" 1 ) 
the bias is modest in HIFLUGCS, groups and poor clus- 
ters appear to be very strongly affected by this effect. Our 
analysis of the sample qualitatively confirms this result. 
While only 12 objects (18%) exhibit a temperature be- 
low 3 keV in the complete sample, almost half (5 out of 
11) of the objects which were rejected based on our flux 
cut in the 0.2rsoo — rsoo annulus fall into this category. 
As a result, an artificial over-representation of cool cores 
among low- luminosity objects might be observed when an- 
alyzing such a sample. Such a dependen ce of the CC frac - 
tion was i ndeed observed in HI FLUGCS ( |Chen et al.|2007 |. 
Recently, Johnson et al. ( 2009 1 claimed to have identified a 
population of merging, NCC galaxy groups, which present 
very shallow surface-brightness profiles, and hence would 
be difficult to detect in the RASS. This population of NCC 
groups could have been missed when selecting the sample, 
thus influencing the observed CC fraction. 

A similar, albeit more modest effect impacts on the 
high-redshift clusters (z > 0.1). At these redshifts, only the 
brightest objects are selected, and because of the cut-off 
in the luminosity function these objects are very rare (see 
Sect. [3|. Our sele ction again qualitatively confirms this re- 
sult (see Sect. 5.3 ). While only three objects in the complete 
sample have a redshift higher than 0.1, two of them (A2204 
and RX J1504) are strong CC objects and are rejected 
when performing the selection based on the flux in the 
0.2r 50 o — fsoo radial range, the rem aining one (A2163) being 
a ver y disturbed merging object ( Markevitch & Vikhlinin 
2001 1. While this result is not statistically significant given 
the very small number of objects in this category, it agrees 
with the predictions of our simulation. This implies that 
the CC bias should be taken properly into account when 
studying the cosmological evolution of the CC fraction. 

In addition, we have shown that our fluxes in the 
0.2r5oo — ?"5oo radial range can be used as a tracer of the 
state of a cluster when comp ared to the core flux integrated 
up to 0.2r 500 (see Sect. [13, since the ratio F annulus / F core 
correlates with the central cooling time. This relation can 
be interpreted as the fact that, while the profiles are very 
different in the cores, they a re essentially self-similar in 



the outer regions of clusters (Neumann & Arnaud 2001 
Leccardi et al. 2010). Indeed, while in the central regions 



non-gravitational effects (radiative cooling, AGN feedback, 
...) are important, in cluster outskirts the profiles are mostly 
determined by gravitational processes. Highlighting the im- 
portance of e xcluding the cor es when analyzing gravita- 
tional effects, |Maughan| ( |2007[ ) observed that the M - T 
relation tightens when these quantities are measured be- 
yond 0.15r5oo- The ratio F annu i us / F core therefore traces the 
deviations from the expectations of the self-similar model, 
and thus can be used to trace the state of a cluster. This 
suggests that for the extraction of quantities depending on 
gravitational processes only, such as the M — T and Lx — T 
relations, the use of a sample selected using fluxes in a well- 
defin ed radial range excluding the cores is important (see 
e.g., |Ota et aL|2006[ ). 



7. Conclusion 

In this paper, we presented an analysis of the effect of the 
different surface-brightness profiles for CC and NCC clus- 
ters on the selection of X-ray flux-limited samples, based 
on both numerical and observational approaches, and esti- 



mated the fraction of CC objects in the nearby Universe. 
Our results can be summarized as follows: 

— Performing realistic simulations of a population of clus- 
ters and of the selection process of the HIFLUGCS sam- 
ple, we estimate that ~ 29% of the strong CC objects 
present in the sample should be removed if one wishes 
to measure the fraction of CC vs NCC objects using this 
sample. 

— Analyzing the populations of simulated CC and NCC 
clusters in the Lx — z plane, we see that CC clusters 
populate a broader range in the diagram. In other terms, 
the flux limit is actually different for CC and NCC ob- 
jects when selecting the clusters based on their observed 
flux only. 

— In all cases, we find that low-luminosity objects (groups, 
poor clusters) are much more affected than more lu- 
minous objects. This effect might explain the lack of 
NCC objects in group samples ext racted from RAS S 
data noted by several authors (e.g., Chen et al. 2007). 
This illustrates the importance of taking into account 
this bias when computing the CC fraction. We also find 
a trend of increasing bias with redshift. 

— From our analysis of the surface-brightness profiles of 
all HIFLUGCS clusters, we propose to select clusters 
according to their flux in a well-defined physical radial 
range excluding the core (0.2r 50 o — ^500) ■ Performing a 
new selection according to the flux in this radial range, 
we exclude 13 objects from the original sample, all of 
which present CC characteristics, and extract a subsam- 
ple of HIFLUGCS as free as possible of the CC bias. 

— Less than half (35 — 37%) of the clusters in our subsam- 
ple exhibit strong CC properties. Relating the observed 
CC fraction with the rate of major entropy injection 
events in the local Universe, we give a rough estimate 
of ~ 0.3 Gyr^ 1 for the rate of heating events per dark- 
matter halo and per Gyr. 

— In addition, we find that the ratio F annu i us / F core be- 
tween the fluxes in the radial range 0.2r 50 o — r 500 and 
— 0.2r5oo strongly correlates with the central cool- 
ing time, and therefore can be used as a CC indicator. 
We expect this indicator to be particularly effective for 
intermediate-redshift objects (z ~ 0.1). 
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Table 4. Basic properties of HIFLUGCS clusters. Column descrip tion: 1: Cluster name ; 2: ins trument used (l=ROSAT; 
2=XMM-Newton; 3=combination of the two); 3: redshift, from |Reiprich fc Bohringcr (2002) and r eferences therein; 4: 
virial temperature in keV, from H10; 5 and 6: rsoo from the scaling relations of jArnaud et ah (2005) in kpc (5), and the 
subtended angle #500 in arcmin (6). 



Cluster 


Instrument 


Redshift 


kT v ir [kev] 


rsoo [kpc] 


$500 [arcmin] 


A85 


3 


0.0556 


6 


1208 


19.2 


A119 


3 


0.044 


5.73 


1187 


23.5 


A133 


3 


0.0569 


3.96 


981 


15.2 


NGC 507 


3 


0.0165 


1.44 


539 


27.5 


A262 


3 


0.0161 


2.44 


729 


38.1 


A400 


3 


0.024 


2.26 


695 


24.6 


A399 


2 


0.0715 


6.7 


1268 


15.9 


A401 


3 


0.0748 


8.51 


1427 


17.2 


A3112 


3 


0.075 


4.73 


1064 


12.8 


Fornax 


1 


0.0046 


1.34 


520 


93.9 


2A 0335 


3 


0.0349 


3.53 


935 


23.1 


Zw III 54 


2 


0.0311 


2.5 


734 


20.2 


A3158 


3 


0.059 


4.99 


1100 


16.5 


A478 


3 


0.09 


7.34 


1316 


13.4 


NGC 1550 


2 


0.0123 


1.34 


519 


35.3 


EXO 0422 


2 


0.039 


2.93 


801 


17.8 


A3266 


3 


0.0594 


9.45 


1514 


22.6 


A496 


3 


0.0328 


4.86 


1098 


28.7 


A3376 


3 


0.0455 


3.8 


966 


18.5 


A3391 


1 


0.0531 


5.77 


1186 


19.6 


A3395S 


3 


0.0498 


4.82 


1086 


19.1 


A576 


2 


0.0381 


4.09 


1005 


22.8 


A754 


3 


0.0528 


11.13 


1648 


27.4 


Hydra A 


3 


0.0538 


3.45 


874 


14.3 


A1060 


1 


0.0114 


3.16 


846 


62.1 


A1367 


1 


0.0216 


3.58 


947 


37.1 


MKW 4 


1 


0.02 


2.01 


651 


27.5 


Zw CI 1215 


2 


0.075 


6.27 


1225 


14.7 


NGC 4636 


3 


0.0037 


0.9 


415 


93.0 


A3526 


1 


0.0103 


3.92 


996 


80.8 


A1644 


2 


0.0474 


5.09 


1117 


20.6 


A1650 


2 


0.0845 


5.81 


1174 


12.7 


A1651 


3 


0.086 


6.34 


1226 


13.0 


Coma 


1 


0.0232 


9.15 


1513 


55.3 


NGC 5044 


1 


0.009 


1.22 


492 


45.6 


A1736 


2 


0.0461 


3.12 


828 


15.7 


A3558 


3 


0.048 


4.95 


1101 


20.0 


A3562 


3 


0.0499 


4.43 


1041 


18.3 


A3571 


3 


0.0397 


7 


1314 


28.6 


A1795 


3 


0.0616 


6.08 


1213 


17.5 


A3581 


2 


0.0214 


1.97 


644 


25.5 


MKW 8 


2 


0.027 


3 


816 


26.0 


RX J 1504 


2 


0.2153 


9.53 


1414 


6.9 


A2029 


3 


0.0767 


8.26 


1405 


16.5 


A2052 


3 


0.0348 


3.35 


866 


21.4 


MKW 3 


3 


0.045 


3.9 


979 


18.9 


A2065 


2 


0.0721 


5.4 


1138 


14.2 


A2063 


3 


0.0354 


3.77 


966 


23.5 


A2142 


1 


0.0899 


8.4 


1408 


14.4 


A2147 


1 


0.0351 


4.07 


1004 


24.6 


A2163 


2 


0.201 


15.91 


1840 


9.5 


A2199 


3 


0.0302 


4.37 


1042 


29.5 


A2204 


2 


0.152 


8.92 


1411 


9.1 


A2244 


1 


0.097 


5.78 


1165 


11.1 


A2256 


3 


0.0601 


7.61 


1358 


20.0 


A2255 


3 


0.08 


5.81 


1176 


13.3 


A3667 


3 


0.056 


6.39 


1247 


19.6 


Sersic 159-03 


3 


0.058 


2.57 


737 


11.2 


A2589 


3 


0.0416 


3.89 


979 


20.4 


A2597 


3 


0.0852 


4.05 


980 


10.5 


A2634 


1 


0.0312 


3.19 


844 


23.2 


A2657 


3 


0.0404 


3.52 


932 


20.0 


A4038 


3 


0.0283 


3.14 


837 


25.3 


A4059 


3 


0.046 


4.22 


1018 


19.3 



14 



Eckert, D. et al.: The Cool-Core Bias in X-ray Galaxy Cluster Samples I 



Table 5. Results of surface-brightness profile fitting. Column description: 1: Cluster name; 2: /?; 3: outer core radius r c \ 
(in kpc); 4: inner core radius r c i (in kpc); 5: ratio between the two beta components at r — 0, R; 6 and 7: estimated 
total fluxes from to r 500 in the 0.5-2.0 keV band, in units of 10~ n ergs cm -2 s _1 , from XMM-Newton /MOS2 (6) and 
ROSAT /PSPC (7); NB: a cross-calibration factor of 0.85 must be applied to the XMM-Newton fluxes (see Appendix[A} 
; 8: minimum x 2 of the fit versus number of degrees of freedom. 



Cluster 


B 






R 




Fx MM 


* / 1 \.y o . i / 


X 2 /d.o.f. 


A85 


0.638±0.009 


192±8 


35±2 


11.2±0.5 


5 


,175±0. 


.072 


4.540±0.038 


349.2/222 


A119 


0.647±0.015 


363±11 






2 


,834±0. 


.076 


2.451±0.042 


309.8/229 


A133 


0.579±0.009 


145±11 


30±1 


23.8±2.7 


1 


.626±0. 


.038 


1.392±0.023 


302.9/220 


NGC 507 


0.458±0.003 


16.4±0.5 






1 


.625±0. 


.036 


1.461±0.034 


784.5/261 


A262 


0.531±0.008 


85±3 


8.0±0.3 


22.6±0.9 


5 


,622±0. 


.187 


4.522±0.115 


428.6/229 


A400 


0.491±0.008 


94±4 






1 


,616±0. 


.070 


1.568±0.031 


339.4/229 


A399 


0.599±0.031 


210±14 






2 


,111±0. 


.085 




101.7/96 


A401 


0.679±0.009 


222±5 






3 


.372±0. 


.048 


3.098±0.034 


313.8/212 


A3112 


0.608±0.006 


88±4 


23±1 


6.6±0.6 


2 


,222±0. 


.023 


1.855±0.021 


263.9/219 


Fornax 


0.685±0.033 


135±4 


12.3±0.9 


8.4±0.4 








11.980±1.062 


501.9/142 


2A 0335 


0.659±0.015 


125±7 


26.0±0.8 


52±7.2 


7 


.000±0. 


.150 


5.655±0.166 


285.4/203 


Zw III 54 


0.703±0.029 


99±5 






1 


.092±0. 


.081 




115.1/115 


A3158 


0.632±0.01 


172±4 






2.41±0.039 


2.196±0.033 


313/193 


A478 


0.662±0.004 


121±2 


31.1±0.9 


3.7±0.1 


3 


,394±0. 


.020 


3.028±0.013 


592.4/229 


NGC 1550 


0.436±0.012 


21±4 


2.8±0.2 


36.1±3.1 


4 


.206±0. 


.501 




127.8/114 


EXO 0422 


0.635±0.016 


71±6 


21±1 


7.5±0.9 


1 


.734±0. 


.067 




148.7/114 


A3266 


0.783±0.013 


413±9 






4 


,061±0. 


.040 


3.298±0.036 


795.5/246 


A496 


0.544±0.006 


96±5 


18±0.2 


17.9±0.9 


6 


,545±0. 


.138 


5.614±0.115 


200.8/144 


A3376 


1.26±0.077 


625±31 






1 


.344±0. 


.069 


1.339±0.044 


137.9/115 


A3391 


0.566±0.015 


163±10 












1.371±0.071 


201.8/146 


A3395S 


0.484±0.023 


173±14 






1 


,465±0. 


.155 


1.311±0.124 


273.2/172 


A576 


0.57±0.013 


107±6 






1 


,911±0. 


.079 




179.7/108 


A754 


0.848±0.015 


406±8 






5 


.488±0. 


.078 


4.368±0.048 


259.8/116 


Hydra A 


0.802±0.009 


173±3 


38.9±0.6 


11.8±0.2 


3 


,177±0. 


.021 


2.741±0.015 


607.9/211 


A1060 


0.667±0.021 


101±7 


27±3 


1.8±0.2 








6.391±0.431 


163.4/144 


A1367 


0.619±0.017 


290±10 












5.077±0.104 


353.1/146 


MKW 4 


0.76±0.076 


186±25 


41±6 


6.2±1 








1.101±0.192 


130.2/112 


Zw CI 1215 


0.652±0.021 


209±9 






1 


.522±0. 


.053 




136.8/116 


NGC 4636 


0.548±0.003 


3.1±0.1 


0.7±0.04 


2.6±0.2 


1 


,571±0. 


.045 


1.361±0.043 


1951. /259 


A3526 


0.633±0.02 


149±5 


20.0±0.9 


16.9±0.8 








16.16±1.145 


275.4/144 


A 1644 


0.448±0.031 


122±23 


17±4 


3.6±1.5 


3 


.043±0. 


.197 




141.8/107 


A1650 


0.866±0.045 


276±20 


83±5 


4.5±0.3 


1 


.494±0. 


.056 




157.4/114 


A1651 


0.637±0.008 


134±4 






1 


,793±0. 


.027 


1.487±0.020 


219.3/159 


Coma 


0.651±0.005 


249±3 












22±0.173 


397/146 


NGC 5044 


0.529±0.002 


8.5±0.1 












4.000±0.048 


657.1/146 


A1736 


0.67 4 


360±4 






1 


,746±0. 


.074 




223.6/107 


A3558 


0.551±0.003 


152±2 






4 


498±0. 


.025 


4.020±0.023 


637.5/216 


A3562 


0.461±0.003 


69±2 






2.12±0.019 


1.734±0.025 


557.5/248 


A3571 


0.639±0.01 


167±8 


58±5 


1.2±0.2 


8 


,227±0. 


.206 


7.417±0.139 


259.5/229 


A1795 


0.724±0.008 


210±5 


54±1 


14±0.5 


4 


,718±0. 


.033 


3.853±0.022 


304.7/199 


A3581 


0.569±0.004 


18.8±0.4 






2 


.083±0. 


.038 




242.5/116 


MKW 8 


0.498±0.051 


94±13 






1 


.362±0. 


.490 




128.3/113 


RX J 1504 


0.759±0.018 


161±12 


53±2 


11.4±1.6 


1 


,459±0. 


.018 




101.4/79 


A2029 


0.616±0.009 


111±17 


43±6 


3.6±1.8 


4 


,967±0. 


.254 


4.159±0.158 


195.3/210 


A2052 


0.749±0.017 


159±4 


31.9±0.8 


15.7±0.5 


3 


.461±0. 


.074 


2.847±0.060 


313.8/229 


MKW 3 


0.628±0.009 


86±5 


27±2 


3.6±0.4 


2 


.486±0. 


.070 


2.014±0.053 


338.1/229 


A2065 


0.755±0.086 


327±50 


93±13 


4.4±0.5 


1 


,716±0. 


.212 




155.2/111 


A2063 


0.734±0.027 


194±13 


54±4 


3.6±0.3 


2 


,597±0. 


.155 


2.330±0.116 


266.2/223 


A2142 


0.757±0.021 


425±28 


140±5 


9±0.9 








3.821±0.033 


148.3/99 


A2147 


0.369±0.01 


61±5 












3.921±0.125 


307.4/190 


A2163 


0.628±0.016 


253±12 






1 


,465±0. 


.029 




89.3/72 


A2199 


0.625±0.004 


99±3 


31±1 


4.6±0.3 


7 


,696±0. 


.084 


6.376±0.049 


393.1/229 


A2204 


0.663±0.013 


134±9 


31±1 


12.7±1.3 


1 


.874±0. 


.037 




145.3/114 


A2244 


0.627±0.011 


103±5 












1.229±0.026 


106.5/86 


A2256 


0.814±0.015 


391±9 






3 


,954±0. 


.069 


3.498±0.036 


527.5/231 


A2255 


0.795±0.018 


424±13 






1 


,264±0. 


.029 


1.054±0.014 


300.6/207 


A3667 


0.529±0.004 


199±3 






5 


,024±0. 


.035 


4.157±0.027 


936.3/261 


Sersic 159-03 


0.733±0.02 


118±13 


40±2 


13.7±2.8 


1 


.718±0. 


.035 


1.444±0.025 


433.7/228 


A2589 


0.625±0.01 


102±4 


21±3 


1.4±0.1 


1 


.824±0. 


.051 


1.583±0.035 


246.1/229 


A2597 


0.708±0.016 


138±11 


41±1 


20.9±2.6 


1 


,545±0. 


.022 


1.256±0.017 


315.8/184 


A2634 


0.626±0.037 


259±21 












1.398±0.086 


158.2/116 


A2657 


0.539±0.005 


77±2 






1 


.717±0. 


.047 


1.611±0.020 


461.3/230 


A4038 


0.53±0.004 


39.7±0.7 






4 


.085±0. 


.056 


3.591±0.045 


286.7/216 


A4059 


0.577±0.006 


71±2 


13±1 


2.2±0.1 


2 


,416±0. 


.038 


1.965±0.031 


259.5/224 



4 See Appendix [b] 
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— ^oo radial range for 37 of the 64 HIFLUGCS clusters, 
for which we performed the analysis by combining both 
instruments. The surface-brightness profiles were fitted si- 
multaneously with the same physical parameters for both 
instruments, while only the normalizations and background 
levels were adjusted individually. F or th e details of the flux 
reconstruc tion procedure, see Sect. |4.2| 

Figure |A.1| shows the comparison between PSPC and 
MOS2 fluxes in the 0.5-2.0 keV band. The left panel shows 
the PSPC fluxes as a function of the MOS2 fluxes, fitted by 
a simple linear relationship. The right panel shows the dis- 
tribution of the ratio between the MOS2 and PSPC fluxes. 
The best-fit relation gives 
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,(A.l) 



Frosat = (0.854±0.004)Fx MM + (2.4±l.l)xl0- 

where the fluxes are expressed in units of ergs cm~ 2 s . 
The scatter of the relation is 5%. 

In conclusion, we can see that a systematic difference in 
absolute calibration is clearly present, XMM-Newton giving 
15% higher fluxes with respect to ROSAT. 



Appendix B: Notes on individual clusters 

— A85: This cluster is categorized as SCC but appears 
to be merging with at least one group South of the core 
(Kempner et al. 2002). This substructure was excised 
when extracting the surface-brightness profile. 



NGC 507: The surface-brightness profile of this nearby 
group (z = 0.0165) is not well represented by either a 
single beta or a double beta model (x 2 = 784.5/261 
d.o.f. for the best fit with a single beta model). 

A 399/ A401: These two cluster s appear to be con- 
nected ( |Sakelliou k Ponman|2004) . The ROSAT /PSPC 
pointed observation was pointed on the centre of A401. 
The surface-brightness profile for A401 was therefore 
extracted in a sector excluding A399. 

Fornax: This nearby poor cluster shows a main peak 
on the BCG NGC 1399 and a secondary peak on 
NGC 1404. NGC 1399 was chosen as the center of 
the cluster and the area surrounding NGC 1404 was 
ignored when extracting the surface-brightness profile. 
Overall, even when using the large ROSAT FOV only 
a small fraction of r 500 is observed, hence the errors 
when extrapolating the fluxes can be large. Moreover, 
the surface-brightness profile of the source is not well 
fitted by any simple model, so the results of the fitting 
procedure might be unstable. We consider this object 
as a very peculiar case. 



Appendix A: XMM-Newton/MOS2 vs 
ROSAT/PSPC cross-calibration 

The present work offers an excellent opportunity to study 
the flux cross-calibration between the instruments used 
(XMM-Newton /MOS2 and ROSAT/PSPC). Indeed, the 
HIFLUGCS sample contains bright objects which are in- 
trinsically persistent over very long timescales, unlike the 
majority of X-ray emitting sources. In this appendix, we 
present a comparison between the unabsorbed fluxes in the 



NGC 1550: This very nearby group (z = 0.0123) 
was observed only by XMM-Newton. As a result, less 
than half of r 500 was observed within the FOV of 
the instrument, so the extrapolated fluxes are rather 
uncertain. 

A3 376: A very disturbed cluster which shows an 
elongated, cometar y shape, and two bow-shock-like 



giant radio relics (Bagchi et al. 2006); no obvious 



center can be defined. Using the emission centroid as 
the center, this cluster shows an anomalously high core 



16 



Eckert, D. et al.: The Cool-Core Bias in X-ray Galaxy Cluster Samples I 




10 



20 30 40 50 60 70 
XMM-Newton/MOS2 flux [ergs cm" 2 s~ 1 ] 



80 



90 



x10 _1 




Fig. A.l. Comparison between XMM-Newton /MOS2 and ROSAT /PSPC unabsorbed fluxes in the 0.5-2.0 keV band. 
Left: PSPC fluxes vs MOS2 fluxes fitted by a simple linear relationship. Right: Ratio between the MOS2 and PSPC 
fluxes. 



radius (r c \ = 625 ± 31 kpc) and an extremely steep 
decline (/? = 1.26 ± 0.08). As a result, the model gives 
only a rough estimate of the surface brightness. 

A 339 5s: This cluster appears to be merging with a 
small er structure in the East (A3395e, Donnelly et al. 
2001). The surface-brightness profile was extracted in 



a sector excluding the secondary structure. 

A 754 '■ The morphology of t his well-known mer ging 

for 



2004 



cluster is very disturbed (see Henry et al.[ 
a detailed analysis of the XMM-Newton data), so no 
obvious center can be defined. We used the emission 
centroid as the center. 



profile. 

A4038: It is not clear whether the surface-brightness 
profile of this cluster should be modeled with a single 
beta or a double beta. A better fit is achieved when 
using a double beta model, but the improvement to the 
fit is not significant, so the single beta model was used. 



NGC 4636: The surface-brightness profile of this 
nearby elliptical in the outskirts of the Virgo cluster is 
not well represented by any simple model (% 2 = 1951.7 
for 259 d.o.f.). As a result, the output parameters are 
uncertain. 



A 1644'- This cluster s hows a double-peak structure 
(Reiprich et al. 2004 1. We extracted the surface- 
brightness profile m a sector avoiding the North 
substructure. 



A1136: This is a difficut case, since no ROSAT/PSPC 
observation exists and the available XMM-Newton 
observation was affected by a high background level. As 
a result, the f3 value could not be constrained, and was 
fixed to the standard value of 0.67 while fitting. We note 
however that our output fl uxes are in good agreemen t 
with the original flux from ( Reiprich & Bohringer 2002 1 . 



MKW 8: This cluster extends well beyond the FOV 
of XMM-Newton, which introduces rather large errors 
when extrapolating the fluxes to r5oo- 

A2634: The ROSAT/PSPC image shows a secondary 
peak North-West of the main cluster. This structure 
was excised for the extraction of the surface-brightness 



